Challenges for scaling cosmologies 
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A cosmological model that aims at solving the coincidence problem should show that dark energy 
and dark matter follow the same scaling solution from some time onward. At the same time, the 
model should contain a sufficiently long matter-dominated epoch that takes place before acceleration 
in order to guarantee a decelerated epoch and structure formation. So a successful cosmological 
model requires the occurrence of a sequence of epochs, namely a radiation era, a matter-dominated 
era and a final accelerated scaling attractor with Q v ~ 0.7. In this paper we derive the generic 
form of a scalar-field Lagrangian that possesses scaling solutions in the case where the coupling 
Q between dark energy and dark matter is a free function of the field ip. We then show, rather 
surprisingly, that the aforementioned sequence of epochs cannot occur for a vast class of generalized 
coupled scalar field Lagrangians that includes, to our knowledge, all scaling models in the current 
literature. 
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I. INTRODUCTION 

The unexpected discovery of the accelerated expansion 
of the universe opened a Pandora's box of new issues and 
questions. Many of these are related to the nature of dark 
energy and to its role within the particle physics model 
(see Ref. pj for reviews). Other questions arise because of 
the so-called coincidence problem: why two components 
that are completely unrelated and scale with time in a 
different way, namely dark energy and matter, appear 
to have roughly the same energy density just now and 
only now ("now" here means within the last one or two 
c- folding times). 

It is possible that once we know the fundamental na- 
ture of dark energy the problem of coincidence will be au- 
tomatically and naturally explained. On the other hand, 
the reverse could be true as well: understanding the ori- 
gin of the coincidence could shed light on the nature of 
dark energy and its relation to the rest of the world. This 
is the footpath that we intend to pursue in this paper. 

This work rests on a fundamental assumption: a com- 
plete solution of the coincidence problem requires that 
dark energy and matter follow the same evolution with 
time, at least from some time onward. For otherwise it 
is clear that the occurrence of coincidence will always 
depend on the initial conditions of the system: chang- 
ing the ratio dark energy/matter at some initial time 
will always imply a displacement in time of the coinci- 
dence epoch. In other words, we can explain the coin- 
cidence only if we show that it is not a coincidence at 
all, but rather that energy and matter always (or from 
some time onward) shared a similar fraction of the total 
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budget. From a phase-space perspective, explaining the 
coincidence requires showing that our present universe 
has already attained its final attractor solution. A solu- 
tion in which matter and dark energy densities are both 
finite and have a constant ratio is denoted in literature 
a scaling solution [2y3j,|J] or, if stable and accelerated, a 
stationary solution |5j- 

Once one accepts this assumption, then it follows the 
first immediate consequence: if we require that the en- 
ergy density of dark energy is proportional to that of 
matter (i.e., prmj p m = const) and at the same time we 
require that the dark energy equation of state parameter 
is less than —1/3 to get an accelerated expansion, then 
one needs either to assume that matter has the same 
negative equation of state as dark energy or that there 
is an interaction between the two components, so that 
p m does not scale as a -3 . The first possibility is clearly 
to be ruled out because such a modified matter equation 
of state would profoundly affect the growth of pertur- 
bations. In fact, for any good model it is not enough 
to require present acceleration: we also need the uni- 
verse to pass through a decelerated matter dominated 
epoch in the past in order to have a well-behaved epoch 
of structure formation. A successful cosmological model 
should therefore admit for a sequence of epochs: a radi- 
ation era, a sufficiently long matter dominated era and a 
final stable accelerated scaling solution. This paper aims 
at searching for such a "good scaling cosmology" . 

The assumption of a stable accelerated scaling solution 
requires therefore the existence of an interaction between 
dark energy and matter [(| . If dark energy is modeled as 
a scalar field then the interaction with matter has to be a 
scalar force additional to gravity. That is, our model has 
to be a scalar-tensor gravitational theory or, equiv- 
alently, an Einsteinian theory with an explicit coupling 
between matter and field. These models have been stud- 
ied for many times in the past and several important 
properties have been discussed [j| 0, . It is known, for 
instance, that a standard scalar-tensor model with an ex- 
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ponential potential has a stable scaling solution and that 
the scaling solution can be accelerated However, in 
this case, it is possible to show that no matter phase pre- 
cedes the acceleration. In other words, after the radiation 
dominated era the system enters directly the accelerated 
regime. This is in contrast with observations, as it has 
been shown in Ref . [|J . 

The fact that the simplest case does not work is the 
main motivation for us to look further. A simple gener- 
alization of the scalar field Lagrangian is to consider the 
so-called fc-essence Lagrangian such that the La- 
grangian is a function p(X, ip) of the field (p and of the 
kinetic term X = — (l/2)g^d IJi <p d u (p. This form is rela- 
tively simple, being still second-order in the field, and has 
been already investigated for several times. We note that 
this type of Lagrangian also covers a wide variety of dark 
energy models such as quintessence ^3$, tachyons ITi). 
phantoms [l^ and (dilatonic) ghost condensates [l6Lll7j . 
If dp/dX < 0, it has been shown that a phantom behav- 
ior (hide < —1) occurs 0, and that matter feels a 
repulsive scalar force [l^. Moreover, it has been shown 
also that if the system contains a scaling solution then p 
can be cast in the form [rH l2p| 

p(X,<p)=Xg(Y), (1) 

where g(Y) is any function of the argument Y = Xe Xip 
with A being a constant. For instance, g(Y) = 1 — c/Y 
is in fact the standard Lagrangian with an exponential 
potential (p = X — ce~ Xv ). The above form for p was 
shown to be valid for uncoupled dark energy and for the 
case in which the coupling 
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is a constant (here S m is the action for matter and gM is 
the metric determinant). 

In this paper we perform a search of a good scaling 
cosmology in three steps. First, we show that the La- 
grangian extends also to the case of variable coupling 
Q(<p) up to a field redefinition. This is an interesting 
result in itself since it unifies some sporadic results ob- 
tained in different ways in literature (e.g. |2l|). Then, 
assuming as a model for g a polynomial 

g = Y J C n Y-\ (3) 

n 

with both positive and negative integer powers of n, we 
derive the critical points of the system. Finally, we show 
that within this class of models there is no way to obtain 
a sequence of a matter phase followed by a stable accel- 
erated scaling solution. When a kinetic scaling matter- 
dominated era exists, this stage is generally followed by 
a scalar-field dominated attractor (f2 v = 1) instead of 
an accelerated scaling attractor. Although our proof of 
absence of two scaling regimes does not extend to any 
possible g(Y), we believe that it seriously undermines 
the real possibility of realizing such an ideal cosmology. 



This negative result opens the challenge: is there any 
case in which a successful sequence can be realized? In 
the final section we will comment on the possibility of ob- 
taining a good scaling cosmology with a fractional power 
Lagrangian g = cq — cY~ u , in which < u < 1. Never- 
theless, we leave a more complete study of this, perhaps 
very exotic, case to future work. 

Beside being a way to approach the problem of coin- 
cidence, a scaling cosmology also provides us with a use- 
ful alternative to standard dark energy scenarios. The 
behavior of the background cosmology and of its linear 
perturbations is in fact radically different in scaling cos- 
mology with respect to most other models. Let us just 
mention three basic differences (see e.g., Refs. 0, H2> ) ■ 
First, in a scaling cosmology the acceleration could start 
at any epoch in the past. Second, the perturbations can 
keep growing even during the accelerated regime. Third, 
the amount of dark energy does not necessarily become 
negligible at high redshifts. All three features radically 
distinguish scaling cosmologies from usual dark energy 
models which only focus on dark energy itself and not 
on its relation with matter. As such, scaling cosmologies 
may serve as a useful testing ground for observations. 

Before passing to the actual calculations, we should 
spend a note on the local gravity constraints on scalar 
forces. In principle, the coupling we introduce is severely 
constrained by local gravity experiments on scalar-tensor 
theories. However, these can be escaped at least in three 
ways. First, by designing a potential with a large mass 
and, consequently, a short interaction range [241 ]. Sec- 
ond, by building a model that happens to satisfy the 
constraints now, but not in the past. Third, by assum- 
ing that the baryons are actually uncoupled to the scalar 
field . The first two solutions change the potential and 
affect the global evolution and therefore will not in gen- 
eral satisfy the requisite for a cosmology that solves the 
coincidence problem. The third case on the contrary can 
be implemented without affecting the potential of the 
scalar field. 

In general, a component of uncoupled baryons can 
dominate in the past, even if their abundance now is very 
small J5(- In this case a matter phase does exist, but it is 
a baryonic matter epoch instead of a dark matter epoch. 
This raises many problems on its own. For instance, the 
baryonic perturbations are almost erased on small scales 
due to the coupling to radiation and therefore, without 
the support from dark matter, would hardly grow to the 
observed amplitude; moreover, the baryonic era would 
finish early in the past, at redshifts quite larger than 1 
and the subsequent accelerated regime would be too long 
to be in fair agreement with both the supernovae exper- 
iments (although here the discrepancy is marginal, see 
Ref. [12]) and with the integrated Sachs- Wolfe effect 
In any case, if a standard (dark) matter phase exists, the 
baryons would never dominate, as we will argue later on. 
Hence the search for a "good scaling cosmology" can sim- 
ply neglect the small baryon component and this is what 
will be done in the present work. Finally, since the two 
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matter components have a different coupling, one has to 
choose a physical frame in which baryons arc conserved 
(otherwise particle masses will be time- varying) but dark 
matter is not. We will work therefore in this frame, which 
is the so-called Einstein frame. 



II. THE GENERAL LAGRANGIAN FOR 
SCALING SOLUTIONS WITH ARBITRARY 
COUPLING 

We shall derive here the general Lagrangian admitting 
scaling solutions in the case where the coupling between 
dark energy and dark matter depends upon the field ip. 
This is the generalization of the works [171, Let us 
start by considering the following action, written in the 
Einstein frame: 



S = d xy/-g M 



M 2 



S m [(P,1pi,9iiu] 



(4) 

where Mp is the reduced Planck mass, R is the Ricci 
scalar, ip is a scalar field, X = — (l/2)g tlv d il tp d v ip and 
ipi are the various matter fields. Notice that we allow 
for an arbitrary coupling between the matter fields and 
the scalar field ip. As mentioned in the Introduction, in 
order to cope with current observations, we assume that 
ip couple only to dark matter. We will also suppose that 
the dark matter component dominates over any other 
baryonic form of matter. 

We are interested in cosmological scaling solutions in a 
spatially flat Friedmann-Robertson- Walker (FRW) back- 
ground metric with a scale factor a(t): 



ds 2 



-dt 2 + a 2 (t)dx 2 



Fricdmann equation in Einstein gravity is given by 
3H 2 = Mp 2 p T , 



(5) 



(0) 



where Mp 2 — 8nG with G being the gravitational con- 
stant, and pt is the total energy density of the universe. 
In what follows we shall set Mp = 1. 

Our focus will be on solutions with constant equation 
of state parameter w v = p(X, ip) / p v in the scaling regime 
and in which the universe is filled only by two compo- 
nents: a barotropic fluid (such that w m = p m /Pm) and 
the scalar field ip. Rewriting the Klein-Gordon equation 
for the field p (in the above metric) in terms of its energy 
density, p v = 2Xdp/dX — p, one gets ^3 



dN 



+ 3(1 + w^pp 



~QPr 



dip 
''dN 



(7) 



where N = In a and Q is defined by Eq. • 

If one starts from scalar-tensor theories Q or a mass- 
varying neutrino scenario [25|. we have 



dp ¥ 



+ 3(1 + w v )p ip = -Q(l - 3w m )p, 



dip 
'dN 



(8) 



This case reduces to Eq. if one relates the coupling Q 
to the coupling Q as Q{\ — 3w m ) = Q. In what follows 
we shall derive the condition for the existence of scaling 
solutions by using Eq. J7J|. Note that the energy density 
of a barotropic fluid satisfies 

+ 3(1 + W m )p m = QPmJ^- (9) 



We shall define the fractional densities of p v and p„ 



as 



_ p v 



3H 2 



«2m — 



3H 2 ' 



(10) 



These satisfy il v + Q m = 1 from Eq. 1(5} . Scaling solu- 
tions are characterized by the condition p v / p m — const, 
in which case fl v is a constant. Using these relations, 
together with Eqs. l[7ll and 10, and following the proce- 
dure of Refs. [13, 123 but dropping the assumption of a 
constant coupling Q(f), one finds the relations 



d In p v d In p v d In p n 



dN 



dN 



dN 



-3(1 + uw), (11) 



and 



dip 
dN 



377 = 7f77_r( m m. — w v ) oc 



1 



Q(v>) v ' Q(v) 

where we introduced an effective quantity: 



(12) 



(13) 



From these equations and the definition of X, one ar- 
rives at 



2X = H 
and thus 



2 / 

dN 



dlnX 
~dN~ 



-3(1 + w cS ) 



d\nQ 
' dN 



(14) 



(15) 



Making use of Eqs. (fTT|> . (fT^|> and l(T5|) we arrive at 
the following generalized "master equation" for the La- 
grangian p(X, ip) : 



2 dQ(ip) 



XQ 2 df 



dhip 1 d In p 



d\nX XQ dip 



= 1, 



where 



l + w m - £l v (w m - ro y ) 



(16) 



(17) 



Equation i|16|) reduces to the one found in Ref. [l7j when 
Q{<p) is constant. 

Solving Eq. (fT7))l one gets: 



p{X,p) = XQ 2 (v)g(xQ 2 {f)> 



(18) 
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where g is an arbitrary function and 



III. PHASE SPACE EQUATIONS 



ip(tp) 



(19) 



See Appendix A for the derivation of Eq. 1)18(1 . In the 
case of constant coupling both Q 2 terms in Eq. I|18|) can 
be absorbed in the definition of g, so our solution reduces 
to that in Ref. [l7^ . In a nutshell, what Eq. (|18|l means 
is that any Lagrangian that allows scaling solutions with 
constant w v can always be cast in the above form by a 
convenient field redefinition. The standard kinetic case 
corresponds therefore to p = XQ 2 — e~ x ^ . Another ex- 
ample is provided by a coupling Q — l/tp as in Ref. H3. 
By using Eq. O we find that the e A ^ term in Eq. JTHJ is 
given by e A ^ — <p x . In this case the Lagrangian (|18fl be- 
comes p — Xg(Xtp x ~ 2 )/p 2 . When A = 2 this simplifies 
to p = g(X)/ip 2 , where g(X) = Xg(X) is an arbitrary 
function of X. This form of p corresponds to the choice 
given in Ref. plf. 

Now let us make the following field redefinition: p — > 
ip(<p), with ip(p) as defined in 119|) . This in turn implies 
X — > = X Q 2 (<p), and Eq. {TBI becomes: 



p(X,p,ip) — X^ g (X$ e 



(20) 



which is the same functional form found in the constant 
coupling scenario [l7|. At the same time the relation 
between S m and the coupling becomes 



1 = 



-1 



5S„ 



(21) 



We have thus shown that the case of a constant coupling 
(Q = 1) is the most general one. In other words, if one is 
interested in scaling solutions, one can always work with 
a Lagrangian in the above form, no matter what kind of 
coupling one has in mind. 

In order to follow the current notation in the litera- 
ture (e.g. [H 13) we will use, instead of ip, the field 
<p(il>) defined by <p(ip) = ip/Q, where Q is a constant. 
In this case one can absorb the Q term that appears in 
the exponential in the argument of g into the definition 
of A. Hence, in what follows, we shall always consider 
the Lagrangian density (dropping the bars on Q and (p) 



p = Xg(Xe^), 



where now A is given by 



A = Q- 



1 



Q v (w m - w v ) 



(22) 



(23) 



It is important at this stage to realize that the system is 
invariant under a simultaneous change of sign of Q and 
A. We can therefore without loss of generality consider 
only the case A > 0. 

In the Appendix B we generalize the above results for a 
more general cosmological background, in which Eq. © 
is replaced by H 2 cx p^. In subsequent sections, though, 
we shall restrict the analysis to Einstein gravity (q = 1). 



So far we have derived the most general Lagrangian 
that possesses scaling solutions. In addition to scaling so- 
lutions there exist other fixed points for the system l|22|) 
characterized by f2 v = 1. In what follows we shall derive 
the autonomous equations taking into account radiation 
to find the general behavior of the solutions. As we men- 
tioned in the Introduction we are interested in searching 
for a good scaling cosmology, namely a sequence com- 
posed of a radiation epoch, a matter-dominated era and 
an accelerated scaling attractor. 

Many results will be shown to hold for any g(Y). How- 
ever we carry out our search assuming as a reference 
model a polynomial expansion in positive and negative 
integer powers: 



Co 



(24) 



n>0 



n'<0 



where Co, c„ and c n < are constants. Note that if we have 
Co = 1 and all other c n , except cj., being zero, our case 
reduces to that of an ordinary scalar field with an expo- 
nential potential |3j . 

For the Lagrangian density (12211 in the presence of pres- 
sureless dust and radiation we obtain the following equa- 
tions 



3H 2 = X(g + 2g x ) + p m + p rad . 



2H = - 



2X(g + gi) + p m + -prad 



(25) 
(26) 



Cp + 3AH(g + gi )<j> + XX[1 - A{g + 2 9l )] 
+AQ Pm = 0, (27) 



where A = (g + bgi + 2g 2 ) 1 and 

g n = Y n d n g/dY T - 



(28) 



The speed of sound, c s , is related to the quantity A by 

mm 



(29) 



When A 1 = the speed of sound diverges. Hence 
no physically acceptable evolution can cross the border 
A- 1 = 0. 

In order to study the dynamics of the above system 
it is convenient to introduce the following dimcnsionless 
quantities: 
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e --V/2 

V3H 



V3H 



(30) 



Then Y is written as Y — x 2 /y 2 . Equation (|25|) gives the 
following constraint equation 



3H 2 



= 1 - M - z 2 . 



where 



x 2 {g + 2 9l ) 



(31) 
(32) 
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It is important to note that, in principle, Q v could as 
well be negative. From Eq. (12(iH we find 



1 dl I, , „ 

i7div=-2( 3 + 3 ^ 2 + z2 )- 



(33) 



By using Eqs. (J57J, and we obtain the au- 

tonomous equations: 



d.r 
dN 



—x 
2 



1 + gx 2 - 2A(g + g x ) + ^z 2 



V6 
2 



A(Q + X)(g + 2 gi )x 2 -Xx 2 



■QA{z 2 



dA 

dz 

dN 



(3 — \/6\x + 'igx 
(-1 + 3.9a; 2 + z 2 ) 

It is useful to notice the relations 
= X(g + 2g 1 ), w v = 

and also 

2x : 



2 + z 2 



9 + 2.9i 



-f + — + tt ip w ip = gx 



(34) 
(35) 
(36) 

(37) 
(38) 



This means that w v > — I for g + gi > and w v < — 1 
for 9 + 91 < 0. From Eq. i|33|) the effective equation of 
state parameter of the system is given by 



2 H 

zip 



gx 



(39) 



Then, in the absence of radiation (z = 0), one has 
w e g — Cl v W(p [see Eq. (|38f) ]. 



IV. CRITICAL POINTS 

In this section we shall derive the fixed points for the 
above autonomous system in the absence of radiation 
(z = 0). The critical points for z 7^ are irrelevant 
to our study. They are listed in the Appendix C, for the 
sake of completeness. The fixed points for z = are de- 
rived by setting dx/dN = and dy/dN = in Eqs. Q3) 
and i|35[l. From Eq. (|33|) we find two distinct classes of 
solutions, either for 3 — V6Ax + igx 2 = or for y = 0. 
The former case gives both scalar-field dominated and 
scaling solutions [l(|. In fact in this case we have 



V6(l + u^f^) 



2A 



(40) 



and, inserting this into Eq. I|34|l we find two cases: 

• Point A: a scalar-field dominated (SFD) solution 
with 



1 



(41) 



• Point B: a scaling solution with 



w v (Q + A) 



(42) 



Let us remind that by definition a scaling solution corre- 
sponds to a situation in which fl v equals neither 1 nor 0. 

The properties of the points A and B will be discussed 
in subsections A and B, respectively. In subsection C 
we shall discuss the second class of solutions, in which 
y = 0. Since these solutions exist in the limit of a van- 
ishing potential, we denote them as kinetic solutions. 



Point A: Scalar-field dominated solutions 



When ftp = 1, we have the following relations ^(j 



9<Xa) = 



\/6Xxa — 3 



3x 2 A 



9i{Y A ) 



6-V6A. 



XA 



6x 2 A 



(43) 



Specifying the model g(Y) one obtains Ya and the 
fixed point {x A ,y A ) by using Eq. (|43J) and the relation 
Y = x 2 /y 2 . From Eq. J2HI we find 



= -1 



V6A 



-XA 



(44) 



A general Lagrangian could have in principle many dif- 
ferent classes of the point A. The number of such critical 
points is known by solving Eq. (|43|l . 

The stability of fixed points can be analyzed by consid- 
ering linear perturbations around them. This was carried 
out in Ref. [l(| for a general g(Y) for positive values of Q 
and A. The eigenvalues of the matrix for perturbations 
are given by 



M+ = -3 + V&(Q + X)x A 



Xx A . (45) 



The fixed point is a stable node if /x+ < and /x_ < 0. 
Allowing negative values of Q as well, the fixed point A 
is stable when 



x A < 



XA <f 
V6 



2(Q+A) 



2(Q+A) ^ JjA ^ A ' 



if Q > -A/2, 

if - A < Q < -A/2. 

if Q < -A. 



(46) 



For negative xa, which corresponds to a phantom 
equation of state (w v < —1) from Eq. 144(1 . the first 
two conditions in Eq. (|46l) are automatically satisfied. 
Hence when Q > — X the phantom fixed point is always 
classically stable. On the other hand the stability of non- 
phantom fixed points (xa > 0) depends upon the values 
of Q and A. Since xa is given by xa = A/[v6(<? + 31)], 
the stability condition l|46|) for non-phantom fixed points 
is expressed as 



.9 + .91 > X(Q - 
.9 + .91 > A 2 /6. 



A)/3, 



if Q > 
if Q< 



-A/2. 
-A/2. 



(47) 
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B. Point B: Scaling solutions 

The scaling solution satisfies the relation (|42ll . Then 
Eq. I0OJ gives 

V6 



XB = WTxy (48) 

We also obtain the following relations valid for all g ^(J : 

2Q(Q + \) 



9{Y B 



\l(p - - 



Q 



Q + X ' 

Q(Q + A) 
Q(Q + X) + 3(g + 9l ) 
Q(Q + X) + 3(g + gi ) 
(Q + A) 2 



(49) 
(50) 

(51) 
(52) 



Again, once the function g(Y) is specified, one obtains Yb 
and db — \xb\/VYb as a function of Q, A. The condition 
for an accelerated expansion corresponds to w c s < — 1/3 
and this gives us 



> A/2 



Q < —X. 



(53) 



which are again independent of the form of g(Y). Note 
that the latter case corresponds to the effective phantom 
(w c s < —1) as we see from Eq. (|5U1) . 

The eigenvalues of the matrix for perturbations around 
the fixed point B are given by |10( 



1 ± Vi - 6 



where 



6 = 
6 = 



3(2Q + A) 
' 4(Q + A) ' 

(1 - n y )(Q + A) 3 [» y (Q + A) + Q] 
3(2Q + A) 2 



(54) 



(55) 



A. (56) 



This point is stable if £i < and £2 > 0. We find that 
negative £1 corresponds to Q > —A/2 or Q < —A. Hence 
when the condition for an acceleration (|53|l is satisfied, 
£1 is automatically negative. In what follows we shall 
consider a realistic situation in which the acceleration 
condition (|53|l is imposed. Then the point B is stable 
when 



Q 



<%<l 

+ A 



and A > . 



(57) 



The second condition is satisfied if we avoid the ultra- 
violet instability of quantum fluctuations 01 > which is 
the case for a non-phantom scalar field. From Eq. 1)521) 
the condition — Q/{Q + X) < Q, v corresponds to —2Q(Q + 
A) < 3(5 + <7i). This is automatically fulfilled for a non- 
phantom fixed point (g + g\ > 0) under the condition 



(|53ll . The most crucial condition for the stability of the 
point B is fl v < 1, i.e., 



g + gi< A(Q + A)/3. 



(58) 



For a non-phantom fixed point this is not satisfied 
if Q < —A but can be satisfied if Q > A/2. Hence 
when Q > A/2 there exist stable, accelerated, and non- 
phantom fixed points B provided that Q v < 1 (whose 
condition is actually required to get a viable scaling so- 
lution). Note that when Q > A/2 the stability condi- 
tion given in Eq. (|4T|> has an opposite equality to that 
in Eq. I|58fl . Hence the stability of the points A and B is 
divided by the border g + gi = X(Q + A) /3, which means 
that the final attractor is either the point A or B de- 
pending on the values of Q and A. When a non-phantom 
scaling solution with positive Q exists in the region 157f) . 
it is the only stable attractor point for any g(Y) , so 
that scaling solutions have the crucial property of being 
global attractors. 

As a last remark, we note that a general form of g could 
in principle exhibit several scaling solutions, all of them 
with the same xb, g(Ys) and w e s, but with different i/b- 



C. Points C and D: Kinetic solutions 

Now we study the second class of solutions of Eq. 
i.e., the case of y = 0. These points exist only if g = 
g{x 2 /y 2 ) is non-singular, i.e., only if one can expand g in 
positive powers of y 2 /x 2 , 



9 = co + ^2> 



n>0 



2 \ n 



(59) 



in which case one has 

g n (y^0) = (n>0). (60) 
In this case Eq. (|34|) is simply given by 

§ = i(3c 0I+ V^)(^-i)=0. (61) 

For Co = this equation gives no real solutions. For 
co 7^ we get the following fixed points: 

• Point C: a t/3-matter dominated era or </?MDE (see 
Ref. 0) 



(xc,yc) 



V6Q 
3c n 



,0 



(62) 



An exception to this rule exists in the case of the fractional 
power-law Lagrangian 1691 with < u < 1, in which a phan- 
tom attractor may also coexist. 
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In this case Eqs. (|^|) and J37J) give 

" v = ~Z , W e g = -5 — , tU<p = 1 • (63) 

3c 3c 

Hence f2 v and Co always have the same sign. When 
c > the solutions is decelerated, and the require- 
ment of the condition Q v < 1 gives 




We note that the cpMDE also corresponds to the 
same class of scaling solutions of point B. In fact 
setting g(Y G ) = c = -2Q(Q + A)/3 in Eq. gSJl 
and eliminating A in Eqs. (|50|l - l|52|) . we obtain the 
results in Eq. JBSJ. Also, we remark that for cq > 0, 
w c s > and therefore dark matter density dilutes 
as a _3 ( 1+u '<=ff) j i.e. faster than baryons. This en- 
sures that baryons did not dominate in the past. 
Since during acceleration baryons dilute faster than 
dark energy, we can safely assume that baryons 
never contributed a large portion of cosmic energy. 
For Co < this is not necessarily true but these 
cases will be ruled out for other reasons presented 
later. 

• Point D: pure kinetic solutions 

(^,^ = (±1/^,0), (65) 

which exists only for positive Co- In this case we 
have that matter is absent and 

Q v = 1 , w cS = 1 , w v = 1 . (66) 

Let us now consider linear perturbations 6x and 8y about 
the generic kinetic fixed point (x,y) — (xh,0). In this 
case the 2x2 matrix M. for perturbations 3] is diagonal 
and its eigenvalues are given by 

3 9 

H\ = -- + -c xl + VGQxk , (67) 

3 ( \/6 \ 

Hi = - ( 1 + c x 2 k - —Xx k j . (68) 

Hence in the case of the ipMDE solution the eigenval- 
ues are fix = (Q 2 / c o) — 3/2 and 112 = 3/2 + Q{Q + A)/co- 
When c > 0, \i\ is negative under the condition 
whereas /i2 > for the values of Q satisfying Eq. (|53|l . 
This shows that the (pMDE corresponds to a saddle point 
for all the relevant cases when Co > 0. When Co is nega- 
tive, it can be a stable point if Q(Q + A) > 3|co|/2. 

In the case of pure kinetic solutions (which exist only 
for Co > 0) one has fii = 3 ± \/&/cq Q and [ii = 
3 =p y/6/^\/2 for x k = ±1/^. Thus, for Q > 0, 
in both cases at least one of the eigenvalues is posi- 
tive, which means that the solutions are either unsta- 
ble nodes or saddle points depending on the values of 
Q or A. When Q < 0, the point x k — Xj^fcQ is stable 
when Q < — -y/3co/2 and A > \/6~co, whereas the point 
Xk = —1/y/co is an unstable node. 



D. Summary of fixed points 

In Table[I] we summarize the property of fixed points 
for the Lagrangian density 1)20(1 . The scalar-field domi- 
nated fixed point A and the scaling solution B exist for 
any form of g(Y) as long as they satisfy the condition of 
existence given in the Tabled Both fixed points can be 
used for late-time acceleration, since the effective equa- 
tion of state w e ff can be smaller than —1/3 depending 
upon the values of Q and A. For a non-phantom case the 
final attractor is either A or B depending on the values 
of Q and A. The scaling solution B is a global attractor 
provided that the condition l|57|) is satisfied. 

The existence of the kinetic fixed points C and D de- 
pends on the form of the scalar-field Lagrangian. They 
appear when g is expanded into positive powers of y 2 fx 1 , 
i.e., Eq. H59|) . An ordinary scalar field with an expo- 
nential potential (g = 1 — c/Y) belongs to this class, 
while for instance a dilatonic ghost condensate model 01 
(g = — 1 + cY) does not. The fixed point C corresponds 
to a saddle point for Co > with f2„ = w c g — 2Q 2 /3cq. 
Hence one can have a temporary scalar-field matter dom- 
inated era (tpMDE) in the presence of the coupling Q. We 
note that Q v < when Co is negative. The fixed point D 
appears only for positive Co and corresponds to fl v = 1 
with no acceleration (w c ff = 1)- Hence this is neither vi- 
able for the matter-dominated era nor for the dark energy 
dominated era and it will not be considered further. 



V. CAN WE HAVE TWO SCALING REGIMES? 

As we anticipated in the Introduction, we search now 
for the occurrence of a two-stage cosmology: a decel- 
erated matter epoch and an accelerated scaling regime. 
This amounts to searching for two distinct fixed points 
for the same set of parameters {Q, A}. Clearly the mat- 
ter point has to be a saddle point in order to give way 
to the final accelerated stable attractor. Since in general 
during the matter epoch there will be a non-negligible 
contribution of the scalar field, this point is, in general, 
a scaling point. Therefore we search for two subsequent 
scaling regimes. It is in principle possible to obtain an 
approximate matter epoch without an associated fixed 
point but this would require a fine tuning of the initial 
condition, so we exclude this possibility here. 

As we have shown in the previous section there are 
two possibilities which lead to an accelerated expansion 
at late times-using either the scalar-field dominated fixed 
point A or the scaling solution B. The <pMDE fixed point 
C appears prior to the accelerated epoch for the models 
given by Eq. (|59|l . For an ordinary scalar field with an 
exponential potential it was found that the (pMDE is fol- 
lowed by the attractor point A (or by point B but in 
this case without acceleration). In this case the present 
universe (Q v ~ 0.7) would be finally dominated by the 
energy density of the scalar field (f^ = 1). Conversely, 
if the present accelerated universe corresponds to a seal- 
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Point 


X 


y 


Existence 


Stability 




w cff 


A 


XA 


(x 2 A /Y A ) 1/2 


ja^O and 17^, = 1 


Stable node under 
conditions 14611 


1 


-1 + sJ%\x A 


B 




(x 2 B /Y B ) 1/2 


3(g + gi) < (Q + A)A 


Stable node for 
-q+a - < 1 and -4 > 


Q(Q+A)+3(g+ S i) 


Q 


2(Q+X) 


(Q+A) 2 


Q+A 


C 


V6Q 
3c 





IQI < V(3co/2) or co < 


Saddle point for Co > 
Stable node for Co < and 
Q(Q + A) > 3|c |/2 


2Q 2 
3c 


2Q 2 
3c 


D 







c > 


Unstable node or saddle for Q > 
Stable node for Q < ~\/3co/2 


1 


1 



Table I: The properties of critical points for the Lagrangian density 1221 in the presence of a pressureless dust (w m = 0). 
Specifying the form of g(Y), xa and Ya are determined by solving Eq. 1431 whereas Yb is known by Eq. 1491 , The kinetic fixed 
points C and D exist when g is given by Eq. 1591 . 



ing attractor B, it was shown that the matter dominated 
epoch, if any, is not sufficiently long to form large-scale 
structure. This is associated with the fact that we re- 
quire a large coupling Q to obtain an accelerated scaling 
attractor, but in this case the solution quickly approaches 
the attractor after the end of a radiation era since there 
is no saddle matter point. Hence one can not have two 
scaling regimes (the decelerated point C and the acceler- 
ated point B) at the same time for the standard scalar 
field with an exponential potential. 

In this section we investigate whether two scaling so- 
lutions can be realized for the general Lagrangian i|2(J|) 
with g given by Eq. I|24|l . We scan all the parameter 
space {Q , A} in search of a successful scaling cosmology. 
Since the cosmological dynamics is different depending 
on the sign of cq, we shall consider three cases (i) cq > 0, 
(ii) Co < and (iii) cq = separately. We shall also 
look into the alternative, fractional power law Lagrangian 
given by 

g(Y) = c Q - cY- u , (69) 
where u, as opposed to n, is not limited to integer values. 



A. Case of Co > 

The function g given in Eq. (|24[1 is composed by pos- 
itive and negative powers of Y. We shall first show that 
the case of positive powers of Y in Eq. (|24Jl is not cosmo- 
logically viable and then proceed to the case of negative 
values of n. 



1 . Positive powers of Y 
Let us first consider the function g given by 

g = c + Y. c « F " n ■ ( 7 °) 

n<0 



In this case all the critical points with y = disappear 
because of the singularity. Then, the only possibilities 
giving rise to a matter-dominated phase corresponds to 
either x — or g — gi — (see Eqs. and J2HJ)), for 
which indeed w c ff = and Q v — 0. For x — > 0, however, 
one has g — * Cq and <7i,<?2 — * from Eq. (|70|l . Then it 
is immediate to see that dx/dN (x — * 0) = —V6Q/(2co) 
from Eq. (|34|) . so we do not have a fixed point unless of 
course Q = (in which case the scaling solution B is not 
accelerated, as can be seen from Eq. (jSJl). If g = g\ = 0, 
the situation is the same and again we only have fixed 
points when Q = 0. Thus we do not have a successful 
cosmological scenario for the function given by Eq. I|7U|I . 
Note also that for the model g = cq —cY~ u with negative 
u the scalar-field energy fraction fl v for the point B has 
to be negative from Eq. H75JI when the condition (|53|l for 
an acceleration is imposed. 



2. Negative powers of Y 

Since we have seen that all positive powers of Y in the 
general polynomial form of g are discarded, let us then 
focus on polynomials with negative power of Y , i.e., 

g = c + CnY~ n ■ (71) 

n>0 

We will prove that when c n =/= for n > 1 it is im- 
possible to have two viable scaling regimes which satisfy 
observational constraints. 

From Eqs. ^ and gSJl we see that the tpMDE de- 
celerated solution C and the accelerated point B have 
always opposite signs of x. In fact, requiring the acceler- 
ation at B, we have either Q > A/2 > or Q < —A < 0. 
In the former case one has xb > and xc < 0, whereas 
in the latter case xb < and xc > 0. However, the 
function g given in Eq. l(7T|l is singular at x — (except 
for power-laws u < 1, see below), which implies that the 
sequence of the solution from C to B is prevented. In 



-1 -0.5 0.5 1 
x 



-1 -0.5 0.5 1 
x 



Figure 1: Phase space for the model l|69^ with u = 2, 
c = co = l,A = 4 and Q = 0.7 together with the fixed 
points A, B, C and D. Here and in the following figures the 
gray area represents the region where Q. v > 1. The dotted 
line corresponds to the singularity given by 1731 at which the 
speed of sound diverges. 



Figure 2: Phase space for the model (1691 with u = 1.1, c = 
co = 1, A = 2 and Q = 2 together with the fixed points A, 
B and D (point C lies in the Q. v > 1 region). The dotted 
line corresponds to the singularity given by 1731 at which the 
speed of sound diverges. 



what follows we will provide a more detailed analysis for 
the possibility of getting two scaling regimes. 

Let us first consider a single power-law function of g(Y) 
given in Eq. I|ti9|) . In the limit x — > with a nonzero value 
of y (which can be y <C 1), the gx 2 y term on the r.h.s. 
of Eq. 1|35[) exhibits a divergence for u > 1 together with 
a divergence of the gx 3 term on the r.h.s. of Eq. I|34[l for 
u > 3/2. In fact, for «/lwe have that 



responds to 



dy/dN 



dx/dN 



(72) 



x^0 



Hence when u ^ 1 the solutions cannot pass the line 
given by x = 0. Since the signs of xb and xc are always 
different, it is inevitable to hit this singularity for u > 1 if 
the solutions move from the i^MDE point C to the scaling 
solution B. This shows that a sequence of solutions from 
C to B is forbidden because of the singularity at x = 0. 

In Fig.^we plot a phase space for the model (|fi9[) with 
u = 2, c = c = 1, A = 4 and Q = 0.7 together with the 
fixed points of the system. The phase space is character- 
ized by another singularity in Eq. I|34|) . associated with 
the divergence of the speed of sound. This appears when 
the quantity, A^ 1 — cq ~ c(u — l)(2u — 1)Y~ U , becomes 
equal to zero, i.e., 



V = ± 



co 



c(u-l)(2u- 1) 



l/2u 



(73) 



For positive c, it exists for u > 1 or < u < 1/2 but 
disappears for 1 /2 < u < 1 . When c < the converse 
is true. 

We note that for the model l|69f) the fixed point B cor- 



xb 



V6 



2(Q + A)' 



VB 



/2Q(Q + A) + 3c \ 1/2M 



3c 



xb 



(74) 

When the condition (|53|l for acceleration is satisfied, we 
require c > for the existence of the point B. Then in 
what follows, we shall only consider the case of positive 
c. When u > 3/2 the point B does not satisfy the con- 

/ \ 1/2" 

dition A > 0, i.e., y < ( ^njp^Tf ) M- This can 
be checked in Fig. 2] in which the scaling solution B ex- 
ists in the region A < 0. When 1 < u < 3/2, it is 
possible to obtain positive values of A. However we still 

t B) 

need another condition: < 1, which gives a more 

severe constraint. Unless u is close to 1, it is not easy for 
the critical point B to fulfill the conditions A > and 

i B) 

fl v < 1. One example satisfying these conditions is 

u = 1.1, Q = 2, A = 2 with Co = ci = 1, as plotted 

in Fig. [21 In this case, however, the <^MDE point exists 
(c) 

in the region > 1. More importantly the trajecto- 
ries can not move from the point C to B because of the 
singularities at x — and also at A^ 1 = 0. 

The above discussion shows that when u > 1 one can 
not realize two scaling regimes. On the other hand, the 
u = 1 case (an ordinary scalar field with an exponential 
potential) is free from both singularities at x — and 
A -1 = 0. This case however has been already ruled out 
as a successful cosmological model @. The argument 
is as follows; we leave u as a free parameter to clarify 
some interesting aspects of the more general case. The 
relevant quantities for the fixed points B and C are given 
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-1 -0.5 




x 



0.5 



Figure 3: Phase space for the model 1)69^ with it = 1, c = 
Co = 1, A = 1.54 and Q = 1.02 together with the fixed points 
A, B, C and D. 



in Eqs. I|50|) . i|63|) and by the following relation: 



(B) _ (2u-l)Q(Q + X) 
* " (Q + A) 2 



3ucq 



(75) 



Let us impose the observational constraints that dur- 
ing the phase C: fl v < 1 and during the phase B: 

/g\ fB) 

w eS < —1/3 and Q v < 1 (notice that the supernovae 
observed value wde is in reality defined through the stan- 
dard Friedmann equation so it cannot be directly used 
here; we have shown elsewhere j^l that the best value 
for itr? is in fact around —0.6 ±0.1). In reality ob- 
servations require quite more stringent constraints than 
this. For instance, supernovae observations constrain 



f£ B) = 0.7 ± 0.2 and 



(B) 
^cff 



< —0.6 ± 0.1. Moreover, too 
much amount of dark energy during the phase C leads to 

a weak growth of perturbations and serious conflicts with 

(c) 

the CMB, so a conservative limit would be f2^ < 0.2 
(see e.g., Ref. |6|). 

The condition for the acceleration of the scaling so- 
lution B {w^' < -1/3) requires cither the condition 
Q > A/2 > or Q < -A < 0. In the latter case it 

(B) 

is easy to find that f2^, ; becomes larger than 1 for u > 1. 
In the former case the condition 3Q > Q + A gives 



v 3Q 2 



2m- 1 



(76) 



while the condition \Q\ < \J3cq/2 for the existence of the 
(ysMDE implies 



8u-3 
9 



(77) 



This shows explicitly that for any u > 3/2 the exis- 
tence/acceleration of the point B is in contradiction with 



the existence of the point C. When u = 1 it is possi- 
ble to have two scaling solutions C and B, but we do 

fB) fB) 

not get values smaller than ur ff = —0.4, O^, = 0.9 

(c) 

and fl v ' = 0.7. This certainly excludes the u = 1 case 
from the range of viable cosmological models. The val- 
ues of the parameters that match this limit are A = 1.54 
and Q = 1.02. The phase space plot in this case is pre- 
sented in Fig. |31 from which it is clear that the ^MDE 
fixed point C is indeed followed by the scaling solution 
B without singularities, although such solutions are not 
cosmologically viable. When 1 < u < 3/2 one can not 
satisfy the observational constraints either, in addition to 
the impossibility of reaching the point B from C. In Fig. [21 
we plot a phase space for the model l|69|) with u = 1.1, 
A = 2 and Q = 2 together with the fixed points of the 
system. 

We have thus shown that it is not possible to obtain 
two ideal scaling regimes for the function (|69|l with u > 1. 
Now we extend this proof to a polynomial form of g(Y) 
with negative powers. The problem as we have seen is 
mainly associated with the fact that two scaling solutions 
B and C are separated by singularities at x = and 
A^ 1 — 0. The latter can disappear by considering the 
sum of the powers given in Eq. (|71|l with the adjustment 
of the coefficients c n . However if the polynomial includes 
any power n larger than 1, this leads to a singularity 
at x = even when the singularity at A^ 1 = is not 
present. Hence if the function g possesses at least one 
term whose power n is larger than 1, the polynomials 
(|71|l are excluded as an ideal scaling cosmology. This 
completes our proof of the impossibility of obtaining two 
scaling solutions in the case of negative powers of Y and 
positive Co- 

Let us conclude this subsection with a brief discus- 
sion of the case u < 1 in Eq. I|69|) . When u < 1 the 
line x = is no longer singular; however, Eq. I|72() still 
holds and the phase space is again separated into positive 
and negative abscissa subspaces. Moreover, the singular- 
ity at A^ 1 = disappears for 1/2 < u < 1 (remember 
we are only interested in c > 0). In this case it is also 
possible to have another nearly matter-dominated phase. 
This corresponds to a situation in which one takes a limit 
x — » with a nonzero but small y in Eqs. I|34() and 13511 . 
Then we can have a matter-dominated era in the region 
x > followed by the scaling solution B with xb > 
(when Q is positive). This situation is similar also in 
the case 0<u<l/2;as long as the system is in the 
region x > and A > initially, the solutions reach 
the scaling attractor B without any singularity. It is in- 
teresting to remark that when u < 1 and c > 0, every 
class A points with xai > is accompanied by a second 
point A: a phantom attractor (xa2 < 0). Hence, it can 
happen that the ipMDE be followed by a point A with- 
out a singularity even for Q > (xc < 0). The u < 1 
"fractional Lagrangians" are then promising but clearly 
for these models to work there are several other obser- 
vational and theoretical issues that should be considered 
and we leave them to a future work. 
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B. Case of Co < 

We shall next consider the case of negative cq. The 
positive powers of Y given in Eq. (|70[) are excluded as vi- 
able cosmological scenarios by arguments similar to those 
presented in the previous subsection. Then let us focus 
on the negative powers of Y given in Eq. J7TJ>- In this 
case one has the t^MDE solution l|62(l with a negative VL V 
[see Eq. 

In addition to the fact that this may be unphysical, we 
are also faced with another problem to obtain two scaling 
regimes C and B. Since the t^MDE satisfies the condition 
(|60|1 one has A^ 1 = g + 5gi + 2g2 = Co, which means 
that A is negative. On the other hand in order to get a 
stable scaling solution B with < 57 v < 1, we require 
A positive. Then, to reach the point B from the point 
C, one needs to cross either the singularity at A^ 1 = 
(which is not allowed) or go through A = (x = 0). 
The latter can only be accomplished in the alternative 
model H69f) with u = 1. Therefore one can not realize a 
good scaling cosmology when Co is negative. 

C. Case of cq — 

The case cq = is also easy to dispose of. In fact in 
this case there are no kinetic solutions and therefore no 
matter eras with y = [see Eqs. (16211 and l|65l) ]. One 
can have a matter era also for x — > 00 with n > 1 or 
for x — > with fractional powers less than one. However 
in both cases A is singular and therefore there are no 
fixed points. Finally, if both y and x go to zero so that 
Y = const., then one can verify that dx/dN does not 
vanish and therefore the point (0, 0) is not a solution. 

This completes our proof. Although the discussion has 
been rather long and technical, the conclusion is straight- 
forward: we have shown that no cosmologically viable 
scaling solutions exist for the general class of integer poly- 
nomial field Lagrangians with variable coupling. 



of fixed points: (A) scalar-field dominated points 
with dp = 1, (B) scaling solution with Q v = 
-Q/[Wtp(Q + A)], (C) a <^MDE solution, and 
(D) pure kinetic solutions. Points of the first two 
classes may exist for any scaling Lagrangian and 
can lead to an accelerated expansion. The acceler- 
ated scaling attractor B, when it exists in the region 
— Q/(Q + A) < fip < 1, is the only global attrac- 
tor apart from the case in which another phantom 
attractor A is present. The points C and D appear 
when the function g can be expanded in the form 
159fl . The yMDE solution C is another scaling so- 
lution (always decelerated in the cases of interest). 

• We have addressed the possibility of finding a se- 
quence of matter and scaling acceleration and found 
that this is impossible for any scaling Lagrangian 
which can be approximated as a polynomial with 
both positive and negative integer powers of its ar- 
gument Y. This is essentially due to the fact that 
a scaling Lagrangian is always singular either along 
the a;-axis or the y-axis of the phase space, thereby 
either preventing the matter-dominated era or iso- 
lating the region with a viable matter era from the 
region where the scaling acceleration occurs. 

It is rather remarkable that the sequence of two scaling 
regimes cannot be realized for such a vast class of scalar- 
field Lagrangians (although, to be fair, we did not in- 
vestigate thoroughly the consequences of Eq. H24[) having 
infinite terms). This underlines how difficult it is to solve 
the problem of coincidence: although cosmological scal- 
ing solutions have been studied for over a decade now, no 
successful case has been identified and this paper shows 
that even a large generalization of the models does not 
help. The search for a good scaling cosmology is not over 
yet, though. In fact we have also shown that a possible 
exception exists in the < u < 1 sector of the function 
g(Y) given in Eq. (|69[) . A detailed investigation of this 
type of fractional Lagrangian is underway. 



VI. CONCLUSIONS 
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Appendix A - Derivation details 

In order to solve Eq. (|16fl . we first rewrite it using as 
a field A times Eq. {Tjft. That is, d/dtp = \Q(tp) d/dip: 



dhip 
dlnX 



2 dQMV)) 



Q(ip) d^ 



dhip 



= 1. 



(78) 
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We then decompose p(X, ip) into 

lnp(X,^) = -^ + ln/(X,V), 
thus arriving at 



background l|86|) . The definition of fJ^ and f2 m are mod- 
ified as 



d In/ 



1 



a In/ 

dtp 



Q(<p(1>)) 

This last equation can be solved by Fourier analysis 
1 



(79) 



0. (80) 



" (#//3 g ) 2 /* ' 



(i///3 9 )2/9 ' 



(87) 



which satisfies + il 



1 from Eq. (EH - 
While Eq. (JTTJ holds for g ^ 1 as well, Eq. (JTIJl is 

subject to the change: 



\nf{X^) 



2X = H 2 



<\nX 



F(u,ip)diJ. (81) 



d</9 
diV 



p«(X,ip) 



Then we obtain the following master equation 



Equation i|80|) then becomes (with = Qoip, where we 
use o to denote a composite function) 



q 



1 



2 dQftp) 



qXQ 2 dip 



d In p 1 9 In p 



dlnX AQ <9v? 



1, (89) 



iujF 



which has as solution 



2 dQ^ 

Qip dtp 



dF 



(82) 



lnF = 



2 dQ^(z) 



Qtp(z) dz 



dz + B(co), (83) 



where A is defined in Eq. I|17|) . The integration of this 
equation gives 

p(X^, J,) = (XQ 2 ( V )) 1/q g {XQ 2 &) e"^) , (90) 

where ip is defined in Eq. (|19fl . For constant Q this re- 
produces the result obtained in Ref. [2(j. 



where B(lu) is an arbitrary function. Undoing the Fourier 
transformation, we get 



In/ = y B(uj) exp 

2 dQ^(z) 
Q^(z) dz 



InX 



dz 



(84) 



-0 + 2 In Q^(tp) + const. 



Finally, recalling l|79|) and noting that the composite 
function Qj, o ijj — Q : we arrive at 



p( x , v) = e 



X e 



Q 2 M 



(85) 



Appendix C - Fixed points for z ^ 

We shall derive here the fixed points for z / 0. In 
this case one has z 2 = 1 — 3<?a; 2 from Eq. i|36|) . Then 
from Eq. I|39|) the effective equation of state always corre- 
sponds to w f[ = 1/3, which means that the scale factor 
evolves as a cx t 1 / 2 . Hence we can not use the fixed 
points with z =/= to get a matter dominated era or 
an accelerated expansion. By substituting the relation 
z 2 = 1 — 3gx 2 for Eq. (|35ll . we find the following two 
cases: (i) y = and (ii) x = 4/(\/6A). 

The case (i) is similar to the kinetic solutions discussed 
in Sec. IV. Then by considering the function g(Y) given 
in Eq. one gets dx/dN = —x(l + y/6Qx) — 0, which 
gives x — or x = — 1/V6Q- Thus, for y — we 
have two fixed points: (a) (x,y,z) = (0,0,1) and (b) 

(x,y,z) = (-1/V&Q, 0,0 - co/2Q 2 ). The point (a) 
corresponds to a standard radiation dominated era with 
fta = 0, whereas for the point (b) there is an energy 
fraction of the scalar field given by f2 v = cq/6Q 2 . 
Appendix B - More general cosmological background In the case (ii) we have X = 4/v / 6A and z 2 = l-8g/X 2 , 

while y is only determined by the specific form of g and 
could as well be zero. From Eq. (|34|l we obtain the 
relation A(X + &Q)(g — gi) = 0, which leads to three 
different fixed points: (c) g = gx, (d) A = — 4Q, and 
(e) A = 0. Solving g — g\ by specifying the form 
of g, we obtain the value y = y c (^ 0), i.e., the fixed 



from which Eq. i|18|) follows immediately by a redefinition 
of the arbitrary function g and by setting V'now - * V'oid/A. 
Notice that the arbitrariness of the function B is ab- 
sorbed into that of g. 



For completeness we also derive the scaling Lagrangian 
in an effective FRW equation which is given by 



#2 _ R 2 n q 



(86) 



where f3 q and q are constants. General Relativ- 
ity, Randall-Sundrum braneworlds j2||, Gauss-Bonnet 
braneworlds [27j and Cardassian Cosmology [2^| corre- 
spond to q = 1, q = 2, q = 2/3 and q = 1/3, respectively. 
The equations Q) and 10 are unchanged even for the 



point (c) (x,y,z) = (4/\/6A, y c , y/1 - Sg/X 2 ). The spe- 
cial case where A = — AQ (i.e., Q < 0) is an interest- 
ing one because yd and Zd are not specified even after 
the form of g is given. In fact, given the form of g, 
the critical curve (d) is found by solving the relation 
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z d = J I - 8g(8/(3A 2 y2)). Finally, the (e) point de- But then z is either imaginary or infinite, so this fixed 

, V, , i . . c ■. 11 point is never realized, 

mands mat y is zero, otherwise since x is timte, g would 

be finite and A = (g + 5gi + 2g-i)~ x would not vanish. 
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